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SUMMARY  - 

Cosmos  1009  rocket  was  launched  on  19  May  1978  into  an  orbit  with  initial 
perigee  height  150  km  and  apogee  1100  km;  its  lifetime  was  only  17  days.  The 
orbit  has  been  determined  daily  during  the  final  14  days  of  its  life,  using  the 
RAE  orbit  refinement  program  PR0P6,  with  about  1100  observations  supplied  by 
NORAD.  An  average  accuracy  of  about  60  m,  radial  and  cross-track,  was  achieved. 

The  orbits  were  analysed  to  reveal  three  features  of  the  upper  atmosphere 
at  heights  between  125  and  175  km.  From  the  decrease  in  perigee  height,  five 
values  of  density  scale  height,  accurate  to  ^4!?,  were  obtained.  The  first  three 
were  within  lOZ  of  those  from  CIRA  1972;  the  fourth,  after  a  magnetic  storm,  was 
higher  than  expected;  the  fifth  gave  evidence  of  the  decrease  in  drag  coefficient 
at  heights  below  130  km. 

Atmospheric  oblateness  produced  a  change  of  4°  in  perigee  ^sition  during 
the  last  four  days  of  the  life.  Analysis  showed  that  the  ellipticity  of  the 
upper  atmosphere  was  approximately  equal  to  that  of  the  Earth,  f  ,  for  the  first 
two  of  the  four  days,  and  about  in  the  last  two. 

The  orbital  inclination  decreased  during  the  14  days  by  about  50  times  its 
standard  deviation,  and  the  observed  variation  was  analysed  to  determine  zonal 
winds  at  heights  of  150-160  km  at  latitudes  near  47  north.  The  zonal  wind  was 
very  weak  (0  i  30  m/s)  for  23-28  May  at  local  times  near  03  h;  and  90  i  30  m/s 
east-to-west  for  29  May  to  4  June  at  local  times  near  01  h. 
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1  INTRODUCTION 

Cosmos  1009  rocket,  1978-50B,  was  launched  on  19  May  1978  into  a  low-perigee 
orbit,  the  initial  orbital  parameters  being  approximately  as  follows:  inclination 
65.1^;  perigee  height  150  km;  apogee  height  1100  km;  eccentricity  0.07;  and 
period  97.4  minutes.  The  rocket  remained  in  orbit  for  only  17  days,  decaying  on 
5  June.  About  1750  observations  were  kindly  supplied  by  the  North  American  Air 
Defense  Command  (NORAD),  of  which  about  1100  were  finally  used  to  determine  daily 
orbits  during  the  final  14  days  in  orbit,  with  the  aid  of  the  RAE  orbit  refinement 
program^,  PR0P6. 

The  orbit  was  chosen  for  analysis  because  its  high  drag  offers  the  chance 
to  determine  upper-atmosphere  properties  at  heights  of  120-180  km  at  close 
intervals  in  time.  Three  distinct  features  of  the  upper  atmosphere  are  investi¬ 
gated.  The  density  scale  height  is  determined  from  the  changes  in  perigee  height 
(sections  3  and  4) .  The  oblateness  of  the  atmosphere  is  assessed  from  the  varia¬ 
tion  of  the  argument  of  perigee  (section  S).  The  zonal  wind  is  studied  by 
analysing  the  decrease  in  orbital  inclination  (section  6). 

2  THE  ORBIT 

About  1100  observations  have  been  used  to  determine  the  orbits,  from  the 
1750  supplied  by  the  assigned  and  contributing  sensors  of  the  North  American  Air 
Defense  Command  (NORAD)  Space  Detection  and  Tracking  Systems  (SPADATS).  Orbits 
have  been  determined  daily  during  the  last  14  days  of  the  satellite's  17-day  life, 
the  epoch  being  at  00  h  each  day.  Table  1  gives  the  orbital  parameters. 

Table  1  shows  that  the  standard  deviation  in  inclination,  i  ,  ranges  from 
0.0003  to  0.0013®;  the  sd  in  eccentricity,  e  ,  from  4  x  10  ^  to  20  x  10  These 
are  equivalent  to  between  30  and  150  m  in  position  for  both  i  and  e  ,  and  the 
average  for  both  is  60  m.  The  sd  in  right  ascension  of  the  node,  0  ,  is  between 
0.0005  and  0.0013®;  while  for  the  argument  of  perigee,  u  ,  and  mean  anomaly  at 
epoch,  ,  the  sd  varies  in  unison,  the  extremes  being  0.005  and  0.054®.  The 
largest  sd  of  each  orbital  element  applies  to  the  final  orbit  (for  the  18  h  before 
decay).  The  inclination  values,  in  Fig  1,  show  a  strong  decrease  which  is  used 
to  determine  the  atmospheric  winds  (section  6) .  The  variations  in  semi  major  axis 
a  and  in  e  and  0  ,  Figs  2-4,  follow  the  expected  pattern.  The  only  unusual 
feature  is  the  'tum-up'  in  u  near  decay  (Fig  5);  this  is  discussed  in 
section  5. 
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3  VARIATION  IN  PERICTIE  HEiqiT 

Perigee  height  is  evaluated  first,  as  a  prerequisite  to  finding  the  heights 
at  which  scale  height,  atmospheric  oblateness  and  atmospheric  rotation  are  being 
measured.  The  values  of  semi  major  axis,  a  ,  and  eccentricity,  e  ,  given  in 
Table  1 ,  were  used  to  calculate  perigee  height  over  a  spherical  Earth,  h^  ,  from 

h_  -  a(I  -  e)  -  R  ,  (1) 


where  R  is  the  Earth's  mean  equatorial  radius,  6378.14  km.  These  values  of  h^ 
are  plotted  as  crosses  with  a  mean  curve  through  them,  in  Fig  6.  The  perturba¬ 
tions,  Ae  ,  in  e  due  to  zonal-harmonic  and  lunisolar  effects  were  calculated 

2 

from  the  PROD  program  ,  and  used  to  obtain  values  of 


h  +  aAe  , 
P 


which  should  be  free  of  gravitational  perturbations.  These  values  of  Q  are 
plotted,  with  a  mean  cuirve,  in  Fig  7.  The  slope  Q  of  this  curve  is  used  to 
determine  density  scale  height  (section  4). 

The  perigee  height,  y  ,  is  obtained  from  (I)  by  adding  the  change  due  to 

P  Q 

the  local  Earth  radius  at  perigee  latitude  to  give,  for  inclination  65.1  , 

y  ■  h  +17.60  sin^u  .  (3) 

P  P 

The  values  of  yp  ,  after  a  small  adjustment  for  the  dynamical  effects  of  the 
Earth's  oblateness^,  are  plotted  as  circles  in  Fig  6,  with  a  mean  (broken)  curve; 
since  u  always  remains  between  52  and  56^,  the  two  curves  in  Fig  6  are  nearly 
parallel. 

4  DENSITY  SCALE  HEICMT 

The  density  scale  height  H  is  a  measure  of  the  rate  of  decrease  of 
density  p  as  height  y  increases,  and  is  defined  by 


1  dp 
”  7  dy 


If,  as  here,  ae/H  >  3  ,  values  of  H  can  be  calculated  from  Q  ,  the  rate  of 
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change  of  Q  due  to  air  drag  using  the  equation 


1 


.la/, . ,  .iL- 2£i  • 

3Mje  y  *  4ae  e 


sin  i  cos  2a) 


where  e'  is  the  ellipticity  of  the  atmosphere  (taken  equal  to  the  ellipticity 

of  the  Earth,  0.00335)  and  H,  applies  at  a  height  y,  ,  1.5H  above  perigee 

I  Ip 

height  .  The  remaining  variables  are  given  in  Table  1. 

Values  of  Q  are  obtained  in  the  form  &Q/At  ,  where  AQ  is  the  change 
in  Q  over  a  suitable  time  interval  At  .  k  change,  AQ  ,  of  at  least  3  km  is 
required  to  ensure  that  the  error  in  Q  does  not  exceed  5  per  cent  and  this 
determines  the  time  interval  At  selected.  Mean  values  of  the  variables  of 
equation  (4)  are  calculated  for  each  time  interval.  In  this  way,  five  observa¬ 
tional  values  of  H|  were  obtained. 

However,  the  method  is  not  satisfactory  for  the  last  2  days  before  decay, 

a 

because  the  values  of  Mj  and  e  change  greatly  between  one  orbit  and  the  next, 
and  the  allocation  of  a  mean  value  for  use  in  equation  (5)  thus  becomes  subject 
to  error,  of  a  magnitude  difficult  to  estimate.  In  these  circumstances  an 
integrated  version  of  equation  (5)  is  preferable.  From  equation  (2)  of  Ref  4, 

Ha  ^1  /3H|  9  \  «, 

h  "  I  +  ^  ~  -  e’a  sin  i  cos  2u)j  ^  ,  (6) 

where  x  >  ae  .  Taking  a  as  constant  (•>  a)  in  the  small  terms  on  the  right- 
hand  side,  taking  cos  2u  as  constant,  and  integrating  between  x^  and  x^  , 
we  have 


H  X  H 

-  X|)  -  -  Xo)  -  T  1”  -  T  (»1  -  V 


a  sin  i  cos  2a) 


If  Xq  and  X|  are  taken  as  the  values  of  x  at  the  beginning  and  end  of  the 
chosen  time  interval  At  ,  Che  left-hand  side  of  equation  (7)  is  AQ  .  The  In 
term  on  the  right-hand  side  is  the  ouiin  term  and  we  may  write  aj  ■  a^  ■  a  in 
the  smaller  terms,  so  that  equation  (7)  gives: 


- 


\  8a 


sin^i  cos  2a) 
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where  e^  and  e^  should,  strictly,  be  cleared  of  zonal  harmonic  and  lunisolar 
perturbations,  so  that  e^  ■  1  -  (R  +  Qq)/*q  »  I978-50B  this  effect  was 

significant  only  for  the  final  value  of  Hj  . 

It  should  be  emphasized  that  equation  (8)  is  preferable  to  equation  (5) 
whenever  (a)  accurate  values  of  e  are  available,  and  (b)  the  drag  substan¬ 
tially  reduces  e  during  the  chosen  time  interval.  At  .  If  e  is  not  much 
reduced  by  drag,  so  that  ^  ,  it  is  dangerous  to  use  equation  (8),  and 

equation  (5)  is  greatly  to  be  preferred. 

For  the  first  three  of  the  values  of  ,  virtually  identical  results  were 
obtained  from  equations  (5)  and  (8);  but  for  the  last  two  values  there  were 
substantial  differences,  and  the  values  from  equation  (8)  are  to  be  preferred  in 
this  situation  for  the  reasons  given  above. 

The  five  observational  values  of  Hj  are  plotted  in  Fig  8  as  circles  with 
a  vertical  error  bar  and  a  horizontal  line  to  show  the  time  over  which  the  value 
applies.  Also  shown  are  the  daily  values  of  the  geomagnetic  planetary  index  Ap 
(plotted  with  a  delay  of  1  day) ,  the  height  y^  ,  and  values  of  Hj  obtained^ 
from  CIRA  1972  for  the  appropriate  height  and  exospheric  temperature,  .  The 
CIRA  values  are  shown  as  crosses.  The  first  three  values  of  H|  are  all  slightly 
higher  than  the  CIRA  values,  the  mean  difference  being  6  per  cent  and  the  maximun 
10  per  cent.  The  last  two  values  are,  however,  16  and  25  per  cent  higher  than 
indicated  by  CIRAt  and  it  seems  that  the  magnetic  storm  on  2  June  1978,  with 
Ap  =  82,  had  a  much  greater  effect  than  is  predicted  by  CIRA.  The  rms  error  in 
H|  is  calculated  to  be  near  I  km  for  all  values,  as  shown  in  Fig  8. 

The  high  value  of  for  3  June  from  the  orbital  analysis  can  be  attri¬ 

buted  to  the  magnetic  storm,  but  the  high  value  on  4  June  probably  arises  from  a 
quite  different  effect  -  a  change  in  the  aerodynamic  flow  regime  at  perigee  from 
free-molecule  flow  to  transition  flow.  Between  June  4.0  and  June  S.O  the  perigee 
height  decreased  from  135  km,  where  the  mean  free  path  of  the  molecules  is  20  m, 
to  126  km,  where  the  mean  free  path  is  7  m.  The  length  of  Cosmos  1009  rocket  is 
believed  to  be  8  m,  and  free-molecule  flow  applies^  only  when  the  mean  free  path 
is  considerably  greater  than  the  dimensions  of  the  satellite.  So  it  cannot  be 
assuned  that  free-molecule  flow  exists  below  about  135  km.  In  transition  flow  the 
drag  coefficient  may  decrease  to  about  half  its  free-molecule-flow  value^.  On 
the  last  day  in  orbit,  therefore,  the  drag  coefficient  would  be  expected  to 
decrease  over  a  short  arc  of  the  orbit  near  perigee,  and,  since  the  orbital 
theory  is  developed  on  the  assumption  of  a  constant  drag  coefficient,  the  result¬ 
ing  decrease  in  drag  would  be  misinterpreted  by  the  theory  as  a  decrease  in  the 
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density  at  perigee,  which  would  be  modelled  as  an  increase  in  scale  height.  Thus 
the  high  value  of  on  the  last  day  represents  a  'drag  scale  height'  rather 

than  a  density  scale  height,  and  would  be  expected  to  exceed  the  value  from  CIRA, 

The  first  four  values  of  H|  have  also  been  plotted  against  height  y^  as 
crosses,  in  Fig  9,  against  a  background  of  comparative  density  scale  height 
curves  for  various  exospheric  temperatures,  T^  obtained  from  CIRA  1972. 

Adjacent  to  each  cross  is  a  value  of  T^  in  brackets,  calculated  using  appropri¬ 
ate  values  of  solar  10.7  cm  radiation  energy,  and  geomagnetic  planetary 

index,  with  adjustments  for  semi-annual  variations.  During  May  and  June  1978,  the 

—22  —2  —1 

value  of  SjQ  ^  was  near  140  x  I0  W  m  Hz  ,  solar  activity  was  fairly  high. 
5  ASSESSMENT  OF  UPPER-ATMOSPHERE  OBLATENESS 

The  argument  of  perigee  u  usually  changes  steadily  throughout  a  satel¬ 
lite's  life,  the  J2  gravitational  perturbation  being  dominant.  However,  if 
u  is  near  45°  (or  135°,  225°  or  315°),  the  pblateness  of  the  atmosphere  creates 
an  asymmetry  between  the  drag  experienced  when  approaching  perigee  and  the  drag 
when  departing  from  perigee,  and  this  asymmetry  has  the  effect  of  altering  u  . 

The  change  in  u  thus  produced  is  too  small  to  be  accurately  measured,  except 
for  satellites  of  very  high  drag,  and  a  substantial  departure  of  u  from  its 
gravitational  variation  is  usually  only  found  during  the  last  day  or  two  in 
orbit.  The  only  previous  analysis  of  this  effect^  was  for  1970-1 14F,  for  which 
the  perigee  height  was  110  km,  and  the  total  change  in  m  was  less  than  1°.  As 
Fig  5  shows,  the  effect  is  much  larger  for  1978-50B,  the  final  value  of  u  being 
about  4°  higher  than  if  the  normal  decrease  had  continued.  It  should  therefore 
be  possible  to  assess  the  magnitude  of  the  atmospheric  ellipticity  c'  . 

For  an  orbit  with  z  *  ae/H  >  3  ,  the  change  in  u>  due  to  atmospheric 
oblateness  is  given  by  equation  (8)  of  Ref  7  as 


Am 

At 


AO  _  .  2e'  sin  i  sin  2m 

-rr  cos  1 - - - 

“  31^.2 


(9) 


where  T.  is  the  orbital  period  expressed  as  a  fraction  of  a  day  and 
“  2 

K  ■  3e'  sin  i  cos  2m  .  The  change  Ad  in  n  which  appears  in  equation  (9)  is 
that  due  to  atmospheric  rotation,  and  equation  (5)  of  Ref  7  shows  that 
AD  ^'I’A  sin  2m  AT^  »  so  that  the  first  term  on  the  right-hand  side  of  equation  (9) 
is  always  less  than  IZ  of  the  second  term. 
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Because  of  the  difficulty  of  assigning  a  mean  value  of  e  over  the  last 
day  or  two  of  the  satellite's  life,  equation  (9)  cannot  be  used  as  it  stands  but 
needs  to  be  integrated.  Equation  (4.52)  of  Ref  8  indicates  that 


“4i[r 


1  + 


-*■  smaller  terms 


s 


(10) 


where  a  “  /T’^^tTt^  and  tj^  is  the  remaining  lifetime  after  time  t^  .  To 
test  the  accuracy  of  equation  (10)  for  1978-50B,  we  have  to  use  the  integrated 
foirm  of  the  equation. 


~  =  I - :p-  0  "  o)  +  In  a  +  smaller  terms  .  (11) 


Fig  10a  shows  the  observational  values  of  T/Tq  (with  the  initial  value  at  MJD 
43660.0,  approximately  4^  days  before  decay),  and  the  variation  given  by 
equation  (11),  taking  t^  =  4.5  days,  which  gives  the  best  fit.  The  agreement  is 
excellent  in  Fig  10,  so  the  theory  is  validated  for  1978-50B  and  the  further 
result  that 


—  ■  a  +  terms  of  0(0.003)  (12) 

®0 
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can  also  be  utilized.  Fig  10b  shows  that  the  decrease  in  e  is  nearly  linear, 
as  predicted  by  (12),  the  slight  increase  in  slope  towards  the  end  probably  being 
caused  by  the  increase  in  density  after  the  oiagnetic  storm  on  2  June.  On  insert¬ 
ing  T  from  (10)  into  (9)  and  writing  e/e^  -  a  ■  (1  -  t/t^)^  ,  we  have 


u  ■  -0  cos  i  + 


.  2. 

sin  1 


sin  2u 


2*o^l 


(13) 


Integration  of  (13),  assuming  sin  2u  ,  K  and  H  are  constant,  gives 


o>  -  u>-  ■  -  40  cos  i  + 


c*  sin  i  sin  2w 


( ^  *  "Ml 

fi  -  -i-V  -  1 

11 

1  h) 

tJ- 

(14) 


d 
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or,  using  (12), 


“-“o' 


=  -  AO  cos  i  + 


2 

e'  sin  i  sin  2u 


(15) 


To  apply  equation  (15),  the  values  of  ai  are  first  cleared  of  gravitational 
perturbations:  the  resulting  values  are  plotted  in  Fig  11a  as  circles.  Also 
shown  in  Fig  11a  are  curves  of  the  variation  of  ui  given  by  equation  (15)  for 
three  values  of  the  atmospheric  ellipticity  e'  ,  namely  0.5f,  0.75f  and  f  , 
where  f  is  the  Earth's  ellipticity  (0.00335).  None  of  the  three  curves  fits 
all  the  points,  and  the  best  fit  is  obtained  by  taking  c*  =  f  for  the  first 
2  days,  and  c'  =  |f  for  the  last  2  days.  Fig  11b.  This  result  suggests  that 
the  magnetic  storm  on  2  June  1978  had  the  effect  of  reducing  the  atmospheric 
ellipticity  to  half  its  normal  value,  and  that  this  effect  continued  over  the 
two  days.  A  decrease  in  oblateness  implies  an  increase  in  density  at  high 
latitudes,  and  this  is  one  of  the  known  effects  of  the  events  associated  with  a 
magnetic  storm. 

In  view  of  this  result,  the  last  two  values  of  scale  height  Hj  were 
recalculated  with  e*  taken  as  ^f  in  equation  (8),  but  the  change  was  not 
significant. 

6  UPPER-ATMOSPHERE  WINDS 

The  14  observational  values  of  inclination,  given  in  Table  1,  were  cleared 

2 

of  lunisolar  and  geopotential  perturbations  using  the  PROD  program  with  numerical 
integration  at  1-day  intervals.  Perturbations  due  to  Earth  and  ocean  tides  are 
expected  to  be  somewhat  smaller  than  the  standard  deviations,  and  were  ignored. 

The  values  of  inclination  cleared  of  perturbations  are  plotted  in  Figs  12  and  13: 
it  is  apparent  that  the  variation  of  inclination  is  very  well-defined,  the 
dec  ^  ^  being  0.03°,  which  is  50  times  the  average  sd  of  the  values. 

The  variation  of  inclination  is  affected  both  by  the  west-to-east  atmos¬ 
pheric  rotation  rate  A  ,  and  by  the  south-to-north  rotation  rate  y  ,  and  the 
theoretical  variation  of  inclination  was  calculated  for  a  series  of  values  of 
A  ,  with  y  *  0  and  y  ■  -0.1,  using  the  RAE  computer  program  ROTATM  with  daily 
integration  steps  between  23  May  and  5  June.  o 

lO 

Unfortunately,  it  happens  that  an  increase  of  0. 1  in  y  has  nearly  the 
same  effect  as  a  decrease  of  0. I  in  A  ,  so  it  is  not  possible  to  exploit  the 
full  accuracy  of  the  values:  instead  we  have  to  give  pairs  of  values  of  A  and 
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y  which  fit  the  variation.  Assuming  that  y  is  fixed,  it  is  not  possible  to 
fit  the  points  with  a  single  value  of  A  ,  and  a  break  has  to  be  made  at  May  29.0. 
Then,  as  Figs  12  and  13  show,  excellent  fittings  are  possible,  for  y  =  0  ,  with 
A  ■  1.0  up  to  May  29.0,  and  A  =  0.8  thereafter;  or,  for  y  ■  -0.1  ,  with 
A  -  0.9  up  to  May  29.0,  and  A  =  0.7  thereafter.  If  the  value  of  y  were 
accurately  known,  the  error  in  A  due  to  observational  error  would  be  about  0.05. 

The  meridional  winds  at  the  relevant  heights  (140-160  km)  and  latitude 
(near  45°  north)  have  been  studied  in  detail  by  Amayenc^,  vdio  determined  the 
steady  component  and  the  diurnal  and  semi-diurnal  variations  in  years  of  high 
solar  activity  (1971-2)  for  all  four  seasons.  His  results  show  that  in  summer 
the  steady  component  is  small,  probably  about  10  m/s  from  north  to  south;  the 
diurnal  component  has  an  amplitude  which  decreases  from  45  m/s  at  160  km  to 
30  m/s  at  140  km,  with  the  maximum  at  21  h  local  time;  the  semi-diurnal  component 
has  an  amplitude  decreasing  from  35  m/s  at  160  km  to  25  m/s  at  140  km,  with  a 
maximum  at  20  h  local  time  at  160  km,  and  22  h  local  time  at  140  km.  For  1978-50B 
these  results  indicate  an  average  south-to-north  wind  of  17  m/s  for  23-28  May 
(height  160  km),  and  a  north- to-south  wind  of  11  m/s  for  29  May  to  4  June  (hei^t 
150  km).  Since  the  likely  errors  in  these  values  are  probably  about  15  m/s,  both 
these  wind  speeds  can  be  regarded  as  small,  and  it  seems  reasonable  to  assume  that 
y  *  0  ±  0.05  (corresponding  to  ±24  m/s).  Thus  the  fitting  shown  in  Fig  12  is  to 
be  preferred. 

Accepting  the  fitting  of  Fig  12,  we  have  A  »  1.0  for  23-28  May  and 
A  *  0.8  for  29  May  to  4  June,  with  sd  made  up  of  0.05  from  the  assumed  error  in 
y  and  0.05  from  observational  error,  ie  a  total  sd  of  0.07.  The  values  apply 
at  a  height  above  perigee,  ie  160  and  150  km  respectively  (on  taking  y^ 

from  Fig  6,  and  H  from  Fig  8).  In  terms  of  wind  speeds,  the  results  are  as 
shown  in  Table  2. 


Table  2 

Values  of  zonal  wind  speed  obtained 


Date 

Latitude 

Local 

time 

Height 

East-to-west 

wind 

Solar 

S|0.7  “■*“ 

Geomagnetic 

conditions 

1978  May  23-28 

48°N 

03  h 

160  km 

0  ±  30  m/s 

145 

Quiet 

May  2 9- June  4 

47°N 

01  h 

150  km 

90  ±  30  m/s 

135 

Disturbed 
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The  latitudes,  local  times  and  heights  quoted  are  of  course  'central*  values:  the 
wind  is  sampled  over  latitudes  between  about  40°  and  55°N,  local  times  vary  by 
up  to  ±2  h  and  heights  by  up  to  ±5  km  from  the  values  given. 

These  results  need  to  be  set  in  perspective  in  relation  to  other  results. 

At  heights  between  120  and  160  km  the  most  numerous  measurements  of  zonal  winds 

are  those  made  by  measuring  the  motion  of  vapour  trails  released  from  rockets, 

which  reveal  great  variability  with  height'^*'*,  and  a  probable  dependence  on 

most  of  the  other  five  parameters  specified  in  Table  2,  namely,  season,  latitude, 

local  time,  solar  activity  and  geomagnetic  conditions.  It  is  also  likely  that 

the  winds  at  these  heights  are  influenced  by  irregular  disturbances  in  the  lower 
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atmosphere,  such  as  weather  systems,  thunderstorms,  etc 

The  zonal  wind  speeds  found  from  incoherent  radar  scatter  measurements,  in 
conjunction  with  theoretical  models,  extend  down  to  130  km,  but  do  not  adequately 
reflect  the  known  variability  of  the  winds  with  height  at  the  lower  levels. 

Results  by  this  technique  for  mid  latitudes  (39-45°)  in  summer  are  given  in 
Refs  13-15;  and  all  indicate  considerable  changes  in  zonal  winds  as  local  time 
changes  from  00  to  04  h,  with  eastward  flow  at  00  h  and  westward  flow  at  04  h. 

The  time  of  zero  wind  is  somewhat  variable,  but  usually  between  02  h  and  03  h, 
which  is  consistent  with  the  first  line  of  Table  2;  the  second  line  of  Table  2 
does  not  show  the  west-to-east  wind  which  is  'normal*  at  01  h,  no  doubt  because 
of  the  occurrence  of  the  geomagnetic  storm,  which  is  quite  enough  to  disrupt  the 
normal  pattern.  (It  is,  of  course,  possible  that  the  meridional  wind  may  also 
be  disturbed:  if  so,  the  resulting  wind  is  likely  to  be  from  north  to  south,  and 
Fig  13  would  then  apply,  ie  the  zonal  wind  would  still  be  from  east  to  west,  but 
stronger.) 

There  are  no  directly  comparable  results  from  satellite  orbit  analysis. 

The  values  at  heights  near  150  km  obtained  by  Forbes^ ^  showed  the  strong  influence 
of  geomagnetic  disturbances  at  latitudes  of  45-75°  north,  but  were  all  for  local 
times  of  10-13  h.  The  eight  values  of  zonal  wind  found  by  King-Hele^^  from 
1969-108A  at  heights  near  160  km  also  showed  the  effect  of  geomagnetic  disturb¬ 
ance,  and  indicated  great  variability  with  local  time.  The  only  one  of  the  eight 
values  at  a  relevant  local  time  is  for  02-05  h  at  160  km  height,  and  gives 
A  ■  0.8  ;  but  it  is  for  equinox  and  for  a  lower  latitude,  20°  north. 

7  CONCLUSIONS 

The  orbit  of  Cosmos  1009  rocket  has  been  determined  at  daily  epochs  during 
the  final  14  days  of  its  17-day  life,  from  about  1100  observations  supplied  by 
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the  North  American  Air  Defense  Command  NORAD.  Though  the  orbit  was  of  extremely 
high  drag,  with  the  apogee  height  decreasing  by  1000  km  in  17  days,  the  perigee 
height  and  inclination  were  determined  with  an  average  accuracy  equivalent  to 
60  m. 

Three  features  of  the  upper  atmosphere  have  been  studied.  The  first  is  the 
density  scale  height.  Five  values  were  calculated  from  the  rate  of  decrease  of 
perigee  height,  and  were  compared  with  values  obtained  from  the  CIRA  1972 
Reference  Atmosphere^.  The  first  three  observational  values  are  all  within  10% 
of  the  CIRA  values,  their  mean  value  being  6%  higher  than  CIRA.  The  fourth  value 
is  16%  higher  than  CIRA:  this  applies  for  the  day  after  a  magnetic  storm,  which 
seems  to  have  had  a  greater  effect  than  CIRA  predicts.  The  fifth  value  applies 
for  the  last  day  in  orbit  ^en  the  perigee  height  decreased  below  130  km,  and 
the  drag  at  perigee  is  reduced  by  the  change  from  free-molecule  to  transition 
flow.  The  quantity  measured  is  then  a  'drag  scale  height',  which  should  (and 
does)  exceed  the  density  scale  height,  being  25%  higher  than  the  value  indicated 
by  CIRA. 

The  oblateness  of  the  upper  atmosphere  produced  a  change  of  4°  in  the 
argument  of  perigee  during  the  final  4  days  of  the  satellite's  life,  and  this 
perturbation  has  been  analysed  to  obtain  values  of  the  atmospheric  ellipticity 
e'  .  The  best  fit  is  obtained  by  taking  e'  equal  to  the  ellipticity  f  of  the 
Earth  (0.00335)  in  the  first  2  days,  and  e'  =  |f  during  the  last  2  days  (June 
3.0  to  June  5.0).  Evidently  the  magnetic  storm  on  2  June  substantially  reduced 
the  atmospheric  oblateness,  presumably  by  increasing  the  density  at  high 
latitudes. 

The  decrease  in  orbital  inclination  is  very  well  determined  (see  Fig  12) 
and  (if  the  meridional  wind  is  weak,  as  expected)  the  observed  decrease  iiiq;>lies 
very  weak  zonal  winds  (0  ±  30  m/s)  at  160  km  height  for  23-28  May,  and  an  east-to- 
west  wind  of  90  ±  30  m/s  at  150  km  for  Nay  29  to  June  4,  when  a  geomagnetic  storm 
occurred.  Both  results  apply  for  latitude  40-50°N,  and  for  local  times  near  03  h 
and  01  h  respectively.  These  results  are  compared  with  values  obtained  by  other 
techniques  in  section  6. 
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The  orbital  inclination  wW  analysed  to  determine  sonal  winds  at  heights  of 
150-160  km  at  latitudes  near  47)|north.  The  sonal  wind  was  very  weak  (0^30  ai/s) 

and  90  ^  30  m/s  east-to-^st  for  29  to 


for  23-26  May  at  local  times  near  03  h; 
4  June  at  locel  times  near  01  h. 


s _ 


Sr  • 


